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1. Introduction 



A central topic in QCD and, more generally, in gauge field theories is the structure of infrared 
(virtual, soft and collinear) singularities of the perturbative scattering amplitudes. The divergent or 
singular behaviour of scattering amplitudes is described by corresponding factorization formulae 
and it is captured by factors that have a high degree of universality or, equivalently, a minimal 
process dependence (i.e. a minimal dependence on the specific scattering amplitude). 

In the following, we consider the factorization formulae describing the singular behaviour 
of QCD amplitudes in kinematical configurations where two or more external-parton momenta 
become collinear. In the case of two collinear partons at the tree level, the collinear-factorization 
formula for QCD squared amplitudes was first derived in Ref. JO]. The corresponding factorization 
for QCD amplitudes was introduced in Refs. ^ ||]. At the tree level, the multiple collinear limit 
of three, four or more partons has been studied [Q, ||, |6[ ^ |8|] for both amplitudes and squared 
amplitudes. In the case of one-loop QCD amplitudes, collinear factorization was introduced in 
Refs. [pi, 10, O, 12], by explicitly treating the collinear limit of two partons. Explicit, though 



partial, results for the triple collinear limit of one-loop amplitudes were presented in Ref. [13]. The 



two-parton collinear limit of two-loop amplitudes was explicitly computed in Refs. [14, 15]. The 



structure of collinear factorization of higher-loop amplitudes is discussed in Refs. [[lfj, |17j, |18| ]. 

The singular collinear factors are customarily expected to depend only on the momenta and 
quantum numbers (flavour, colour, spin) of the collinear partons, with no dependence on the ex- 
ternal non-collinear partons. This feature of collinear factorization is denoted as strict collinear 
factorization, and it is generally assumed to be valid in the calculation of cross-sections in hadron- 
hadron collisions from the convolution of universal (i.e., process independent) parton distribution 
functions with the hard-scattering cross-section. Strict collinear factorization, however, is vio- 



lated beyond the tree-level for amplitudes in space-like collinear configurations [17]. The violation 



is originated [ |17| , |18[ ] by long-wavelength absorptive contributions (such as those produced by 
Coulomb-Glauber gluons [19, 20, 21, 22, 23, 24, E5|]) that causally disconnect initial-state and 



final-state interactions, thus limiting the factorization features due to colour coherence. Owing 
to the absorptive ('imaginary') origin of the violation of strict factorization, the effect is partly 
canceled at the level of squared amplitudes. Indeed, such a cancellation is complete up to the 
next-to-leading order (NLO). Nonetheless, strict factorization is violated at higher orders. This 



challenges the validity of the factorization theorem of mass (collinear) singularities Q26| , |27| ] and 
related issues in the context of the factorization of transverse-momentum dependent distributions 
[28, 29, 3C], and it can produce logarithmically-enhanced radiative corrections [ ^2| , |l7| ] to hard- 
scattering processes in hadron-hadron collisions. 



2. Generalized collinear factorization at all orders 

A set {pi, . . . ,p, n } of m (m > 2) parton momenta approaches the multiparton collinear limit 
when they become parallel. In this limit all the parton subenergies 



sa = (Pi + Pe) 2 , with i,£eC = {l,...,m} 



(2.1) 
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are of the same order and vanish simultaneously [Q, |5]|. The collinear direction is defined through 
the light-like vector 

177 n^ 

rt-rtm-p*—, (2-2) 

where n^ 1 is an auxiliary light-like vector (n 2 = 0), which parametrizes how the collinear limit is 
approached, and p\ m = p\ H Yp m - The longitudinal-momentum fractions zi are 

n ■ pi n-pj 

Zi = — - = , i£C , (2.3) 

n-P n-p Lm 

and they fulfill the constraint Y!i=\ Zi = 1, with p? — > ZiP^{i £ C) in the collinear limit. According 
to our notation, pf is the outgoing momentum in the scattering amplitude, and the time component 
('energy') p° is positive (negative) for a final-state (initial-state) parton. In the time-like (TL) 
collinear region all the parton subenergies in Eq. ( |Q| ) are positive and 1 > z.i > 0; in all the other 
kinematical configurations, we are dealing with the space-like (SL) collinear region. Therefore, in 
the TL case, all the collinear partons are either final-state partons or initial-state partons. In the SL 
case, at least one collinear parton is in the initial state and, necessarily, one or more partons are in 
the final state. The SL collinear limit is typically encountered by considering initial-state radiation 
in hadron collision processes. 

In the multiparton collinear limit the matrix element ^ of the scattering process with n exter- 
nal partons (n > m) fulfills the all-order generalized factorization formula JTv| ] 

Pm > P> Pm+ 1 ) • • • i Pn ) 1^0 , (2.4) 

where Sp is the all-order splitting matrix, which captures the dominant singular behaviour in the 
multiparton collinear limit, and ^# = ^{P,p m+ \ ,p n ) is the reduced matrix element, which is 
obtained from the original matrix element ^ by replacing the m collinear partons A\,... ,A m with 
a single parton A, which carries the momentum P. The matrix element and the splitting matrix 
satisfy the perturbative (loop) expansion: 



+ Jt^ + Jt {2) + . . . , (2.5) 
Sp = Sp(°) + SpW + Sp( 2 ) + ... , (2.6) 

where the superscripts (k) (k = 0, 1,2, . . . ) refer to the order (number of loops) of the perturbative 
expansion. The splitting matrix is expected to be universal and process independent (strictly fac- 
torized), namely, it should depend on the momenta and quantum numbers (flavour, colour, spin) 



of the external collinear partons only. Nonetheless, according to Eq. ([2.4]), the splitting matrix can 
also acquire a dependence on the non-collinear partons. Strict collinear factorization holds at the 
tree-level (in both the TL and SL regions): 

Sp i0) (p u ...,p m ;P;p m+1 ,...,p„) =Sp i0) (pi,...,p m ;P) , (TL and SL coll. lim.) , (2.7) 

and it also holds (because of colour coherence) in the TL collinear region to all orders: 

SpOi,..., )=Sp{p h ...,p m ;P) , (TL coll. lim.) . (2.8) 
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In the SL region, strict collinear factorization is violated [ |T7| ] at one-loop and higher-loop orders. 
We illustrate the violation of strict collinear factorization by mainly considering the infrared 



(IR) divergent part of the splitting matrix. The IR structure of Sp is not independent [ ]13| , [14) , |31| , |32| , 
33] of the IR structure of the QCD amplitude Using dimensional regularization in d = 4 — 2e 
space-time dimensions, the all-order matrix element M fulfills the IR recursion relation [ 34 , 35 



36, 37, 38, 39] 



Me) \M) + |^ fin -} , 

r(2) 



where the operator Im(£) = ijjf '( £ ) ' ( £ ) + ••• (with perturbative coefficients I^ ; (£)) is IR 
divergent, while the matrix element term ..# fin - = „#( ' + JK^)^ 1 - + ^( 2 ' fin - + . . . is IR finite and 
its first contribution in the perturbative expansion is the complete tree-level matrix element 
in Eq. 



(2.9) 



(*), 



An expression analogous to Eq. ( |2.9[ ) holds for j$ and the coiTesponding IR operator 
Ijj-(e). The all-order splitting matrix also fulfills a recursion relation [[L 



Sp = V(e) Sp fin - V 1 (e)= [l-V(e) V _1 (e)] Sp + V(e) Sp fin - V '(e 



-l 



(2.10) 



where the first term in the perturbative expansion of the IR finite splitting matrix Sp fin is the tree- 



level splitting matrix Sp' ' (Sp 1 



.fin. 



S p (0) +Sp (l)fin 

= 1-I(e), 



+ Sp (2)fin - + ...), and 



V _1 (e) 



l-I(e), 



(2.11) 



where I and I are obtained from the collinear limit of the IR operators Im and Ijj, respectively. Up 
to two loops, the expansion of the coefficient of the first term in the right-hand side of Eq. ( |2.1C ) 



reads l-V(e) V -1 (e 



I 



(i) 



I«(e 



(ce|), where 



!(i) 



I 



(i) 



(2.12) 
(2.13) 



The IR operators in Eqs. ( 2.12| ) and ( 2.13| ) have been calculated in Ref. Jl7| ] starting from the 
known IR structure of scattering amplitudes to two-loop order [34, ^]. In the SL collinear region, 
the operator ij^(e), which describes the IR divergent part of the one-loop splitting matrix Sp' 1 ', 
contains factorization breaking contributions that are proportional to the anti-Hermitian operator 



2n 



in 



£ Tj-Ty ®(-Zi) sign(j ) 



(2.14) 



ieC 
jeNC 



where NC = {m + 1, . . . ,n} denotes the set of non-collinear partons, and is the colour-charge 



matrix of the k-th parton (we are using the general colour space notation of Ref. [40]). The operator 



in Eq. ( J2. 14{ ) embodies colour correlations between collinear and non-collinear partons that are 
produced by the non-Abelian Coulomb phase, and thus they violate strict collinear factorization. 
In the two-parton collinear case, these colour correlations are illustrated in Fig. [l] (left). Note that 
the two-parton one-loop SL splitting matrix is known to all orders in e Jl7| ] and, therefore, the 
result of Ref. [JTvJ] is not limited to the treatment of Coulomb-Glauber gluon effects to leading IR 
accuracy. For three or more collinear partons, the IR finite part of Sp' 1 ' is unknown (the explicit 
calculation for three collinear partons is in progress), but it also contains factorization breaking 



contributions similar to those in Eq. (2.14) 
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Figure 1: Two-parton factorization breaking correlations at one-loop in hadron-hadron collisions (left). 
Three-parton factorization breaking correlations at two-loops, n > 4 QCD partons (right). 



At two loops, new factorization breaking terms appear in the SL region through the opera- 



tor [17] 



v (2;2; 

x mC 



271 



1 



nfabcY, I r, a r / fo r, c 0(-z,-)sign(, !7 )0(-^,) 

ieC j,keNC 



x In 



Sjp s kp 
r (2) 



iO 



(2.15) 



,(2; 2) 



which contributes (to V^' c and) to the IR divergent part of Sp' 2 ^. The operator A^'c' includes both 
Hermitian and anti-Hermitian contributions, and it embodies three-parton correlations involving 
one collinear and two non-collinear partons (see Fig. [I] (right)). The Hermitian part also depends 
on the size of the non-collinear momenta. The anti-Hermitian part, which depends only on the sign 
of the partons subenergies, can be rewritten [|l7|] in terms of correlations between two collinear and 
one non-collinear partons. 

The breaking of strict collinear factorization in the SL region beyond the tree-level is due 
to absoiptive contributions originated from the fact that colour coherence is limited by causality, 
which distinguishes initial-state from final-state interactions. Lepton-hadron DIS is, however, a 
special case since all the non-collinear coloured partons are produced in the final state, and thus 
there are no initial-state interactions between collinear and non-collinear partons. The one-loop 



and two-loop SL multiparton splitting matrices explicitely calculated in Ref. Q17[ ] effectively take a 
strictly-factorized form in DIS. 



3. Squared amplitudes and cross-sections 

The all-order singular behaviour of the squared matrix element \^\ 2 (summed over the colours 
and spins of the external partons) in a generic kinematical configuration of m collinear partons is 



obtained by squaring the generalized factorization formula in Eq. ( |2.4[ ) 

~ (Jt\ P(pi,...,p m ;P;p m+h ...,p n ) 



(3.1) 
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where the matrix P is the square of the all-order splitting matrix Sp , 

P=[S P rSp, (3.2) 

and it satisfies a loop expansion similar to Eq. (j2l]|), with perturbative coefficients PW [k = 
0,1,2,...). 

The splitting matrix in hadron-hadron collision processes breaks strict collinear factorization 
beyond the tree-level. However, owing to their absorptive origin, the factorization breaking effects 
partly cancel for squared amplitudes, and thus a violation of collinear factorization in physical 
observables is not directly implied. Factorization breaking contributions can cancel either in the 
squared splitting matrix P, after taking the expectation value with the reduced matrix element in 



Eq. (3.1), or even among different partonic processes with different number of external partons 
contributing to the same physical cross-section at a given order. 

The tree-level collinear matrix p(°) is obviously strictly factorized also in the SL region. The 



divergent part of P^ 1 ) is strictly factorized because the factorization breaking operator in Eq. (2.14) 



is anti-Hermitian. In the two-parton collinear limit also the complete pW, including the IR fi- 
nite part, is strictly factorized [|17[], in spite of the fact that Sp^ 1 ' violates collinear factorization. 



Similarly, the three-parton factorization breaking correlations appearing in the two-parton splitting 
matrix Sp™ (Fig. |] (right)) survive in P^ 2 ' but only for processes with n > 5 QCD partons [O]. 

~(2) — (2) 

The expectation value of the two-loop operator A P (e) = (e) +h.c, which gives the (IR 
dominant) factorization breaking contribution to the squared splitting matrix p( 2 ) for the two-parton 
collinear limit 1 , onto the reduced matrix element j$ is 

Sp^) 1 4 2) ( £ ) SpW = (3# (0) | (Sp(°)) + 4 2) ( £ ) Sp(°) P# (0) > 

| (Sp (0) ) ' 4 2) ( £ ) S P (0) l^ (0) > + c - c - + hi g her orders • ( 3 - 3 ) 



+ 

~(2) 

Although Ap ' is not vanishing, its lowest-order expectation value (i.e., the first term on the right- 



hand side of Eq. (|3.3|)) vanishes in pure QCD [ |18| ] (i.e., if the lowest-order reduced matrix element 
is produced by tree-level QCD interactions). Note however that this contribution would be 
non-vanishing, for instance, for tree-level quark-quark scattering produced by electroweak inter- 
actions (with CP-violating electroweak couplings and/or finite width of the Z and W bosons), or 
if the tree-level QCD scattering is supplemented with one-loop (pure) QED radiative corrections. 

~(2) 

Therefore, as a matter of principle, it remains true that the operator A p explicitly uncovers 
two-loop QCD effects that lead to violation of strict collinear factorization at the squared amplitude 
level. Moreover, the second term in the right-hand side of Eq. ( |3.3| ) is not vanishing and contributes 
to the SL collinear limit of three-loop QCD squared amplitudes, together with the factorization 
breaking effect produced by P^ 3 ' and highlighted independently in Refs. [17, 18]. In the case 



of scattering amplitudes with n = 5 QCD partons, the colour correlation structure of these two 
factorization breaking contributions at three loops is analogous to the commutator structures that 



'The operator A^? (b), which is computed in Ref. p7t|, is analogous to the multiparton collinear operator A^ (e) 
but it includes also the subdominant i/e poles and some finite terms to all orders in e, which are known in the two-parton 
collinear configuration. 
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were found in the N 4 LO computation of super-leading logarithms in 'gaps-between-jets' cross 



sections [22] 



4. Summary 

Collinear factorization has been generalized to all orders for kinematical configurations where 
two or more external partons become collinear, and explicit results on one-loop and two-loop (and 
three-loop) amplitudes for both the two-parton and multiparton collinear limits have been presented 



Q17|]. In the space-like region strict factorization is violated beyond the tree-level for scattering 
amplitudes. For squared amplitudes the strict-factorization breaking terms partly cancel, but the 
remaining effects still lead to logarithmically-enhanced contributions at high perturbative orders 
and challenge the validity of factorization theorems of mass (collinear) singularities. 
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